Alchemical Free Energy Calculations on Membrane-Associated Proteins.

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

Structure-based virtual screening discovers potent and selective adenosine A1 receptor antagonists.

Development of subtype-selective leads is essential in drug discovery campaigns targeting G protein-coupled receptors (GPCRs). Herein, a structure-based virtual screening approach to rationally design subtype-selective ligands was applied to the A1 and A2A adenosine receptors (A1R and A2AR). Crystal structures of these closely related subtypes revealed a non-conserved subpocket in the binding sites that could be exploited to identify A1R selective ligands. A library of 4.6 million compounds was screened computationally against both receptors using molecular docking and 20 A1R selective ligands were predicted. Of these, seven antagonized the A1R with micromolar activities and several compounds displayed slight selectivity for this subtype. Twenty-seven analogs of two discovered scaffolds were designed, resulting in antagonists with nanomolar potency and up to 76-fold A1R-selectivity. Our results show the potential of structure-based virtual screening to guide discovery and optimization of subtype-selective ligands, which could facilitate the development of safer drugs.

Design of Drug Efficacy Guided by Free Energy Simulations of the ß2 -Adrenoceptor.

G-protein-coupled receptors (GPCRs) play important roles in physiological processes and are modulated by drugs that either activate or block signaling. Rational design of the pharmacological efficacy profiles of GPCR ligands could enable the development of more efficient drugs, but is challenging even if high-resolution receptor structures are available. We performed molecular dynamics simulations of the ß2 adrenergic receptor in active and inactive conformations to assess if binding free energy calculations can predict differences in ligand efficacy for closely related compounds. Previously identified ligands were successfully classified into groups with comparable efficacy profiles based on the calculated shift in ligand affinity upon activation. A series of ligands were then predicted and synthesized, leading to the discovery of partial agonists with nanomolar potencies and novel scaffolds. Our results demonstrate that free energy simulations enable design of ligand efficacy and the same approach can be applied to other GPCR drug targets.

Importance of Binding Site Hydration and Flexibility Revealed When Optimizing a Macrocyclic Inhibitor of the Keap1-Nrf2 Protein-Protein Interaction.

Upregulation of the transcription factor Nrf2 by inhibition of the interaction with its negative regulator Keap1 constitutes an opportunity for the treatment of disease caused by oxidative stress. We report a structurally unique series of nanomolar Keap1 inhibitors obtained from a natural product-derived macrocyclic lead. Initial exploration of the structure-activity relationship of the lead, followed by structure-guided optimization, resulted in a 100-fold improvement in inhibitory potency. The macrocyclic core of the nanomolar inhibitors positions three pharmacophore units for productive interactions with key residues of Keap1, including R415, R483, and Y572. Ligand optimization resulted in the displacement of a coordinated water molecule from the Keap1 binding site and a significantly altered thermodynamic profile. In addition, minor reorganizations of R415 and R483 were accompanied by major differences in affinity between ligands. This study therefore indicates the importance of accounting both for the hydration and flexibility of the Keap1 binding site when designing high-affinity ligands.

Adenosine A2A receptor antagonists: from caffeine to selective non-xanthines.

A long evolution of knowledge of the psychostimulant caffeine led in the 1960s to another purine natural product, adenosine and its A2A receptor. Adenosine is a short-lived autocrine/paracrine mediator that acts pharmacologically at four different adenosine receptors in a manner opposite to the pan-antagonist caffeine and serves as an endogenous allostatic regulator. Although detrimental in the developing brain, caffeine appears to be cerebroprotective in aging. Moderate caffeine consumption in adults, except in pregnancy, may also provide benefit in pain, diabetes, and kidney and liver disorders. Inhibition of A2A receptors is one of caffeine's principal effects and we now understand this interaction at the atomic level. The A2A receptor has become a prototypical example of utilizing high-resolution structures of GPCRs for the rational design of chemically diverse drug molecules. The previous focus on discovery of selective A2A receptor antagonists for neurodegenerative diseases has expanded to include immunotherapy for cancer, and clinical trials have ensued. LINKED ARTICLES: This article is part of a themed issue on Structure Guided Pharmacology of Membrane Proteins (BJP 75th Anniversary). To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v179.14/issuetoc.

Fragment-based design of selective GPCR ligands guided by free energy simulations.

Fragment-based drug discovery relies on successful optimization of weakly binding ligands for affinity and selectivity. Herein, we explored strategies for structure-based evolution of fragments binding to a G protein-coupled receptor. Molecular dynamics simulations combined with rigorous free energy calculations guided synthesis of nanomolar ligands with up to >1000-fold improvements of binding affinity and close to 40-fold subtype selectivity.

Structure-Guided Design of G-Protein-Coupled Receptor Polypharmacology.

Many diseases are polygenic and can only be treated efficiently with drugs that modulate multiple targets. However, rational design of compounds with multi-target profiles is rarely pursued because it is considered too difficult, in particular if the drug must enter the central nervous system. Here, a structure-based strategy to identify dual-target ligands of G-protein-coupled receptors is presented. We use this approach to design compounds that both antagonize the A2A adenosine receptor and activate the D2 dopamine receptor, which have excellent potential as antiparkinson drugs. Atomic resolution models of the receptors guided generation of a chemical library with compounds designed to occupy orthosteric and secondary binding pockets in both targets. Structure-based virtual screens identified ten compounds, of which three had affinity for both targets. One of these scaffolds was optimized to nanomolar dual-target activity and showed the predicted pharmacodynamic effect in a rat model of Parkinsonism.

Publisher Correction: GPCRmd uncovers the dynamics of the 3D-GPCRome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

GPCRmd uncovers the dynamics of the 3D-GPCRome.

G-protein-coupled receptors (GPCRs) are involved in numerous physiological processes and are the most frequent targets of approved drugs. The explosion in the number of new three-dimensional (3D) molecular structures of GPCRs (3D-GPCRome) over the last decade has greatly advanced the mechanistic understanding and drug design opportunities for this protein family. Molecular dynamics (MD) simulations have become a widely established technique for exploring the conformational landscape of proteins at an atomic level. However, the analysis and visualization of MD simulations require efficient storage resources and specialized software. Here we present GPCRmd (http://gpcrmd.org/), an online platform that incorporates web-based visualization capabilities as well as a comprehensive and user-friendly analysis toolbox that allows scientists from different disciplines to visualize, analyze and share GPCR MD data. GPCRmd originates from a community-driven effort to create an open, interactive and standardized database of GPCR MD simulations.

Energy Landscapes Reveal Agonist Control of G Protein-Coupled Receptor Activation via Microswitches.

Agonist binding to G protein-coupled receptors (GPCRs) leads to conformational changes in the transmembrane region that activate cytosolic signaling pathways. Although high-resolution structures of different receptor states are available, atomistic details of allosteric signaling across the membrane remain elusive. We calculated free energy landscapes of ß2 adrenergic receptor activation using atomistic molecular dynamics simulations in an optimized string of swarms framework, which shed new light on how microswitches govern the equilibrium between conformational states. Contraction of the extracellular binding site in the presence of the agonist BI-167107 is obligatorily coupled to conformational changes in a connector motif located in the core of the transmembrane region. The connector is probabilistically coupled to the conformation of the intracellular region. An active connector promotes desolvation of a buried cavity, a twist of the conserved NPxxY motif, and an interaction between two conserved tyrosines in transmembrane helices 5 and 7 (Y-Y motif), which lead to a larger population of active-like states at the G protein binding site. This coupling is augmented by protonation of the strongly conserved Asp792.50. The agonist binding site hence communicates with the intracellular region via a cascade of locally connected microswitches. Characterization of these can be used to understand how ligands stabilize distinct receptor states and contribute to development drugs with specific signaling properties. The developed simulation protocol can likely be transferred to other class A GPCRs.

Ligand design by targeting a binding site water.

Solvent reorganization is a major driving force of protein-ligand association, but the contribution of binding site waters to ligand affinity is poorly understood. We investigated how altered interactions with a water network can influence ligand binding to a receptor. A series of ligands of the A2A adenosine receptor, which either interacted with or displaced an ordered binding site water, were studied experimentally and by molecular dynamics simulations. An analog of the endogenous ligand that was unable to hydrogen bond to the ordered water lost affinity and this activity cliff was captured by molecular dynamics simulations. Two compounds designed to displace the ordered water from the binding site were then synthesized and evaluated experimentally, leading to the discovery of an A2A agonist with nanomolar activity. Calculation of the thermodynamic profiles resulting from introducing substituents that interacted with or displaced the ordered water showed that the gain of binding affinity was enthalpy driven. Detailed analysis of the energetics and binding site hydration networks revealed that the enthalpy change was governed by contributions that are commonly neglected in structure-based drug optimization. In particular, simulations suggested that displacement of water from a binding site to the bulk solvent can lead to large energy contributions. Our findings provide insights into the molecular driving forces of protein-ligand binding and strategies for rational drug design.

FZD5 is a Gαq-coupled receptor that exhibits the functional hallmarks of prototypical GPCRs.

Frizzleds (FZDs) are a group of seven transmembrane-spanning (7TM) receptors that belong to class F of the G protein-coupled receptor (GPCR) superfamily. FZDs bind WNT proteins to stimulate diverse signaling cascades involved in embryonic development, stem cell regulation, and adult tissue homeostasis. Frizzled 5 (FZD5) is one of the most studied class F GPCRs that promote the functional inactivation of the ß-catenin destruction complex in response to WNTs. However, whether FZDs function as prototypical GPCRs has been heavily debated and, in particular, FZD5 has not been shown to activate heterotrimeric G proteins. Here, we show that FZD5 exhibited a conformational change after the addition of WNT-5A, which is reminiscent of class A and class B GPCR activation. In addition, we performed several live-cell imaging and spectrometric-based approaches, such as dual-color fluorescence recovery after photobleaching (dcFRAP) and resonance energy transfer (RET)-based assays that demonstrated that FZD5 activated Gαq and its downstream effectors upon stimulation with WNT-5A. Together, these findings suggest that FZD5 is a 7TM receptor with a bona fide GPCR activation profile and suggest novel targets for drug discovery in WNT-FZD signaling.

Agonist-induced dimer dissociation as a macromolecular step in G protein-coupled receptor signaling.

G protein-coupled receptors (GPCRs) constitute the largest family of cell surface receptors. They can exist and act as dimers, but the requirement of dimers for agonist-induced signal initiation and structural dynamics remains largely unknown. Frizzled 6 (FZD6) is a member of Class F GPCRs, which bind WNT proteins to initiate signaling. Here, we show that FZD6 dimerizes and that the dimer interface of FZD6 is formed by the transmembrane α-helices four and five. Most importantly, we present the agonist-induced dissociation/re-association of a GPCR dimer through the use of live cell imaging techniques. Further analysis of a dimerization-impaired FZD6 mutant indicates that dimer dissociation is an integral part of FZD6 signaling to extracellular signal-regulated kinases1/2. The discovery of agonist-dependent dynamics of dimers as an intrinsic process of receptor activation extends our understanding of Class F and other dimerizing GPCRs, offering novel targets for dimer-interfering small molecules.Frizzled 6 (FZD6) is a G protein-coupled receptor (GPCR) involved in several cellular processes. Here, the authors use live cell imaging and spectroscopy to show that FZD6 forms dimers, whose association is regulated by WNT proteins and that dimer dissociation is crucial for FZD6 signaling.

The tyrosine Y2502.39 in Frizzled 4 defines a conserved motif important for structural integrity of the receptor and recruitment of Disheveled.

Frizzleds (FZDs) are unconventional G protein-coupled receptors, which activate diverse intracellular signaling pathways via the phosphoprotein Disheveled (DVL) and heterotrimeric G proteins. The interaction interplay of FZDs with DVL and G proteins is complex, involves different regions of FZD and the potential dynamics are poorly understood. In the present study, we aimed to characterize the function of a highly conserved tyrosine (Y2502.39) in the intracellular loop 1 (IL1) of human FZD4. We have found Y2502.39 to be crucial for DVL2 interaction and DVL2 translocation to the plasma membrane. Mutant FZD4-Y2502.39F, impaired in DVL2 binding, was defective in both ß-catenin-dependent and ß-catenin-independent WNT signaling induced in Xenopus laevis embryos. The same mutant maintained interaction with the heterotrimeric G proteins Gα12 and Gα13 and was able to mediate WNT-induced G protein dissociation and G protein-dependent YAP/TAZ signaling. We conclude from modeling and dynamics simulation efforts that Y2502.39 is important for the structural integrity of the FZD-DVL, but not for the FZD-G protein interface and hypothesize that the interaction network of Y2502.39 and H3484.46 plays a role in specifying downstream signaling pathways induced by the receptor.

Fragment optimization for GPCRs by molecular dynamics free energy calculations: Probing druggable subpockets of the A 2A adenosine receptor binding site.

Fragment-based lead discovery is becoming an increasingly popular strategy for drug discovery. Fragment screening identifies weakly binding compounds that require optimization to become high-affinity leads. As design of leads from fragments is challenging, reliable computational methods to guide optimization would be invaluable. We evaluated using molecular dynamics simulations and the free energy perturbation method (MD/FEP) in fragment optimization for the A2A adenosine receptor, a pharmaceutically relevant G protein-coupled receptor. Optimization of fragments exploring two binding site subpockets was probed by calculating relative binding affinities for 23 adenine derivatives, resulting in strong agreement with experimental data (R2 = 0.78). The predictive power of MD/FEP was significantly better than that of an empirical scoring function. We also demonstrated the potential of the MD/FEP to assess multiple binding modes and to tailor the thermodynamic profile of ligands during optimization. Finally, MD/FEP was applied prospectively to optimize three nonpurine fragments, and predictions for 12 compounds were evaluated experimentally. The direction of the change in binding affinity was correctly predicted in a majority of the cases, and agreement with experiment could be improved with rigorous parameter derivation. The results suggest that MD/FEP will become a powerful tool in structure-driven optimization of fragments to lead candidates.